Nonequilibrium electronic transport through a polymer chain: Role of solitons 
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Nonequilibrium electronic transport through a polymer chain is investigated by the scattering 
state operator method. The polymer chain is described by an electron-lattice coupling model and 
its two ends are connected with metal electrodes of different chemical potentials. The scattering 
states are shown to be a set of complete eigenstates of electrons in the system at nonequilibrium 
steady state. With the method, we show that the nonequilibrium Peierls transition (NEPT) does not 
survive the lattice relaxation and the soliton-antisoliton pair excitations. Furthermore the electronic 
transport though the chain is shown to be accomplished through the soliton-lattice energy band. 

PACS numbers: 73.40.Sx, 72.10.-d, 71.20.Rv 



I. INTRODUCTION 

Molecular electronic devices are attracting great inter- 
ests because of their tremendous advantages in applica- 
tions. A lot of works are focused on electronic transport 
through a sandwich structure. For example, the measure- 
ment of conductance in a single molecule, such as a single 
conjugated polymer wire coupled with metallic leads, is 
recently developed to a more precise level than before us- 
ing many different experimental methods^— To advance 
electronic devices at the single-molecule scale, we are in 
an increasing need for understanding the detail of charge 
transport through an individual molecule. 

Theoretically, a nonadiabatic dynamical evolution 
method has been used in the studies for the polaron for- 
mation and motion in a polymer chaini^r— In these stud- 
ies, the metal electrodes are treated as finite chains due to 
the limitation of numerical computation and then it's not 
applicable to a long-time steady-state. To get a steady- 
state condition, the electrodes are better to be considered 
as (infinite) reservoirs. When the system is staying in a 
steady state with a constant particle/energy flow, i.e. a 
nonequilibrium steady state (NESS), the quantum sta- 
tistical property becomes quite different from that in the 
equilibrium case. It is the special statistical property 
that plays an essential role in the formation of the flow 
and the behaviors of the polymer chain. 

Non-equilibrium Green's function techniques are used 
widely in the mcsoscopic transport problems^ Using 
these methods, noninteracting problems could be solved 
exactly. For systems with many-body interactions, it 
involves summing a set of real time correlation func- 
tions for the linear, quadratic, cubic, etc, response to 
these interactions. Perturbation theory has succeeded 
in many problems before. Today, with great advances 
in experiments, there are more and more problems for 
which the nonperturbative methods are necessary and 
essential. A different approach called scattering state 
operator method, is focused on the the construction of 
steady-state nonequilibrium ensemble*^ In this method, 
nonequilibrium quantum statistical mechanics for steady 
state is written in a form similar to that in equilibrium 
case. Hcrshficlc— has shown the existence of a statis- 



tical operator Y which accounts for the nonequilibrium 
boundary condition as a part of an effective Hamilto- 
nian. And he proposed that Y be constructed in terms 
of the scattering state operators. Once Y operator is 
constructed, the concepts in the equilibrium case and 
the equilibrium numerical techniques can be applied to 
nonequilibrium steady state only if H is replaced by 
H noncq = H — Y . And the well-established equilibrium 
numerical techniques, such as exact diagonalization, the 
quantum Monte Carlo method, renormalization group, 
variational methods, etc, can bring great improvements 
if applied to nonequilibrium steady state. 

Recently, Han et al. have done much based on the 
scattering state operator method^ They established the 
method by explicitly solving the electronic transport 
through a noninteracting single quantum dot (QD). Then 
the method is applied to a single QD with electron- 
phonon or electron-electron interactions and various 
transport phenomena are reproduced which have been 
studies before using other methods such as nonequilib- 
rium Green's function techniques. This work is to gen- 
eralize a single quantum dot to a quantum chain with 
some sites. We solve the transport through a fixed disor- 
dered chain without electron-phonon interactions explic- 
itly, and then apply the method to a polymer chain that 
is described by the Su-Schricffcr-Hcegcr (SSH) models 
with electron-lattice interactions. The work provides a 
ground for the further generalization to the quantum 
chain with many-body interactions. 

We consider the metal electrodes as large reservoirs 
and adopt the scattering state operator method to inves- 
tigate the nonequilibrium electronic transport through 
a polymer chain coupled with the reservoirs. By the 
scattering operator method, we show that nonequilib- 
rium Peierls transition (NEPT) docs not survives the 
lattice relaxation and the soliton-antisoliton pair exci- 
tations, which is in contrary with the result of a previous 
work based on the continuum version of SSH model with 
the order parameter is assumed to be a constant ^ Fur- 
thermore, we show that a soliton lattice^ is formed in 
the NESS and the electronic transport is accomplished 
through the soliton-lattice energy band. 

The arrangement of this article is as follows. In 
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the following section, we give our model for the 
metal/polymer/metal structure and describe the scatter- 
ing state operator method. Numerial results are pre- 
sented in Sec. III. The conclusions of this article are 
given in Sec. IV and details of calculations are given in 
the appendixes. 



where g ak is the coupling between the a electrode and 
the site that is attached at the electrode, here, we fixed 
the site by cum^ = cuN) a , that is, the first (n = 1) site 
of the chain is coupled with the left electrode while the 
last (n = N) site with the right electrode. 



II. MODEL AND METHOD 

A. Model Hamiltonian 

We consider a sandwich structure which consists of a 
polymer chain and two metal electrodes coupled with the 
chain's two ends. The Hamiltonian of the whole system 
is composed of three terms as follows: 



H = H S + H B + V, 



(1) 



where H§ is the Su-Schrieffcr-Hccgcr (SSH) model, which 
describes the polymer chain of iV-sites, 



Hs — H e \ + Hlatt 

with the electronic part as 



(2) 



N-l 



/ j tn f c n,o- c n+l,CT + c n+l,a c n,aj , (3) 



where t n [= to — a (u„+i — u n )\ is the hopping integral be- 
tween sites n and n + 1 with to being the one in an equi- 
distance lattice, a describing the electron-lattice coupling 
between neighboring sites, and u n the monomer displace- 
ment of site n from its equidistant position due to the 
electron-lattice interaction, and a (c n ,a) is the creation 
(annihilation) operator for electron with spin a at site n, 
and the lattice part of the chain as 



„ N-l , r N 

#latt = — K>+1 - Un) +y 

71 — 1 71 — 1 



(4) 



where K is the elastic constant due to the a bonds of 
the polymer chain and M the mass of a site that consti- 
tutes the chain. The metal electrodes are described by 
reservoirs of free electrons 



2^ e akdl kc 

aka 



d a 



(5) 



where a(= L, R) indicates the left and the right metal 
reservoirs and d} aka (d akcT ) is the creation (annihilation) 
operator for electron with energy t ak and spin a in the a 
electrode. The last term in the whole Hamiltonian is the 
coupling between the electrodes and the polymer chain, 
which gives as 



V = ^2 g a k (dl ka c aa + c aa .d ak( /j , 



(6) 



B. Scattering State Operators Method 

The scattering state operator ip aklT satisfies the 
Lippman-Schwingcr equation as follows: 



i 



aka 



e ak - jCo + irj 



(7) 



where rj is an infinitesimal convergence factor. The Li- 
ouville operator Cq is defined as Co A = [Ho, A]. The 
scattering state \ipaka) includes forward and backward 
scattering on the reservoirs and the molecule chain as 
well as the incoming reservoir state \d a ka)- And the 
Lippman-Schwinger equation can be rewritten into the 
Heisenberg's equation of motion as follows: 



H, iplJ = eofcVL + i<lWlka ~ 4k*)- ( 8 ) 



First, we consider a fixed-lattice chain with a set 
of hopping integrals {t n }. The scattering state opera- 
tor ip aka can be expanded within the one-particle basis 

{dlk'o-^La} aS 



^aka = d lka + E ^ + E T&*4'« 



bW 



By inserting Eq.© into Eq.©, we obtain 



Ibk- 



9bk< 



7Z% 



(9) 



(10) 



(here, we have used 7£ fc = 7" and 7^ = 7^,) and 
equations for {7^} 

(eak + iT L )ll k + t\lf = 9ak5 a ,L; 

*„-i7«-i + £afc7f + tnln k +i = (for n = 2, 3, JV - 1); 

*iV-l7jV-l + ( e ak + iTR)lN = 9a k S a ,R, 

(11) 

where we have defined the density of states (DOS) of the 
reserviors 



\9ak\ 



e - £afc + iri 



(12) 



In the large bandwidth limit, the functions T Q (e) are 
energy-independent constants. From the above equa- 
tions, 7^ fe can be written as 



(13) 



then we can obtain the functions i^°^(e), that are given 
in the appendix A. 
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As done by Han^ we can also show the completeness 
of the scattering state operators tp aka defined above un- 
der the assumption that there exist no isolated energy 
cigenstatcs in the NESS, that gives 



JV 



E ^aka^aka = E d lka d ^a + E C n,<T C «,^ (14) 



ak 



a k 



n=l 



which guarantees the original operators {d^^,^ a } can 
be conversely represented by the scattering state opera- 
tors {V'ifco-}: which gives 



(15) 

Furthermore, it can also be shown that the total Hamilto- 
nian in Eq.([T]) can be represented by the scattering state 
operators as follows: 



H = E e ^i ka ^ak 



(16) 



aka 



Then we have the effective Hamiltonian for the NESSii 

H eS = H - Y = (eafe - Ma) i>l k ^ak*, (17) 



which leads immediately 



{^aka^bk'a') 



Tr{e-P H *«} 

= fa{£ak)8a,b5a,cr'5k,k' , (18) 

where (• • • ) stands for the average with respect to the 
NESS, and f a (e) = 1 / \e^~ ^ <> ' kT + 1) is the Fermi dis- 
tribution function, fj, a being chemical potential of the a 
electrode. 



C. Current at Nonequilibrium Steady State 

Now we can calculate the current flow in the sand- 
wich structure. The current from the left electrode to 
the polymer chain is given as usual, 



(19) 



With respect to ip^ aka , the current flow can be evaluated 
and finally represented as the Landauer-Biittiker form 

/=! J deT(e)[f L (e)-f R (e)], (20) 

where T(e) is transmission rate, which can be determined 
by the coefficients {7^ fe } and it gives 



bk 1 

T{e) = AT L T R \F { 1 R \e)\ 2 . (21) 

With the same process, the current between neighboring 
sites in the polymer chain and that from the chain to the 
right electrode can be obtained. With the Eqs. (| 1011 I"3"|) . it 
is easily proved that all them are the same, which meets 
the requirement of a steady state (See Appendix B). 



D. Lattice Relaxation at Nonequilibrium State 

Now we consider the flexible lattice which we treat 
classically in this work. Since the mass of sites in the 
polymer chain is much heavier than that of electrons, the 
Born-Oppenheimer approximation is applicable. Then 
the force exerting on the polymer sites from the electron- 
lattice interaction could be obtained from the Hcllmann- 
Feynman theorem. Furthermore, we consider the system 
to be at zero temperature, then the damping induced by 
exterior environment should eventually make the lattice 
at rest. At the steady state the force exerting on the 
lattice should be canceled out, which gives 



(22) 



where we have defined the dimcnsionlcss staggered lat- 
tice order parameter y n = (— l)™a(u n +i — u n )/to, 
the dimensionless electron-lattice coupling constant 
A = 2a 2 / (irKto), the bond-charge density p n ,n+i = 



4 



,Cn+i,a + h.c), with 8 being a Lagrangian mul- 



tiplier to guarantee the polymer chain length unchanged 
(£„K+i - Un) = 0). By using Eq.®, (USD, and JTg), 
we have 



Pn,n+l - 2^ — 



deKc 



Fi a >(e)Fi%(e)] f a (e). (23) 



III. NUMERICAL RESULTS 

For the sake of simplicity, we consider a symmetric 
case in which the couplings with the electrodes are the 
same, i.e., Ti = Tr = T, and the chemical potentials 
of the two electrodes are taken as = —ftR = V&/2 
with Vb being the applied bias voltage. Furthermore, we 
take to(= 1) as the unit of energy and the dimensionless 
electron-lattice coupling constant A = 0.2 in accord with 
that for £rans-polyacetylene. 



A. Lattice Configuration at Zero Bias Voltage 

Before we present our results on steady current 
through a polymer chain and the role of solitons in it, the 
effects onto the polymer chain arising from the couplings 
with the reservoirs under equilibrium case ( Vb = ) are 
studied. As is well known, the lattice is dimerized due 
to the Peierls instability for an isolated polymer chain. 
The staggered lattice order parameter y n 0.13 except 
a few bonds near its two ends which are affected by the 
boundary condition. The zigzag of the order parameter is 
due to the finite-size effect under the restriction of chain- 
length unchanged. When the couplings to the reservoirs 
are turned on, the lattice configuration is affected by the 
reservoirs as shown in Fig. [T] It is clear that with weak 
couplings, such as T = 0.26 in Fig. QJa), the lattice con- 
figuration should be very similar to that of an isolated 
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FIG. 1: The staggered lattice order parameter {y n } of a 
30-site chain under different couplings with metal electrodes 
without applied bias voltage. The coupling constants are 
taken as T = 0.26 (a), 0.32 (b), and 0.38 (c). The straight 
line indicates the dimerization magnitude at thermodynamic 
limit. 



FIG. 2: The staggered lattice order parameter {y n } of a 30- 
site chain coupled with metal electrodes with a voltage bias 
V b = 0.48 (a), 0.56 (b), and 0.72 (c). The coupling constant 
T = 0.32. 



chain, but a slight reduction of the order parameters near 
the ends, which makes the bond-length changes smaller 
than that for a free boundary condition. Once the cou- 
pling r is increased to an intermediate value T = 0.32 
in Fig. Hlb) , the bond- length changes at the chain ends 
become almost similar as that of the inside bonds. The 
boundary effect on the bond-length at the chain ends is 
canceled out by the couplings with the metal electrodes. 
When the couplings arc enlarged more , such as T = 0.38, 
y n for the bonds near the two ends arc much less than 
0.13, which is opposite to that of an isolated chain with 
a free boundary condition (Fig. HJc)). In the following, 
we take the intermediate coupling T = 0.32. 



B. Lattice Configuration at Finite Bias Voltage 

First, we see the lattice configuration at a finite bias 
voltage, which is shown in Fig. [2] It can be seen that 
the lattice has a little change under a small bias voltage, 
and then with increasing bias voltage, the lattice contains 
one soliton-antisoliton (SS) pair for Vb = 0.56 and two 



SS pairs for Vb = 0.72. It's expected the number of SS 
pairs will be increased one by one and a soliton lattice 
is formed at thermodynamic limit. To understand the 
formation of SS pairs in the polymer chain, we show the 
corresponding transmission rate T(e) in Fig. [3J When 
Vb is below than 0.5, electrons from the reservoirs can 
hardly flow to the polymer chain, due to the gap exist- 
ing in the energy band of the polymer. {y n } shown in 
Fig. [2ja) is similar to that of an isolated polymer chain 
and Tie) is shown in Fig.G^a). When Vb is above a criti- 
cal value, such as Vb = 0.56, the Fermi surface in the left 
reservoir reaches the n = +1 energy level of the poly- 
mer chain while the Fermi surface in the right is below 
n = —1 energy level (see Fig.rj|. Since the whole system 
is on steady state, the n = +1 level is always occupied 
by electrons from the left reservoir and electrons on the 
n = — 1 level flows to the right persistently. In all the 
polymer chain has no extra electrons since the bias volt- 
age is added symmetrically, which corresponds to that an 
electron on n = —1 level is excited to the n = +1 level, 
as shown in Fig. |4jb). In order to minimize the total 
energy of the whole system, the polymer chain will be 
relaxed and a SS pair is formed (Fig.Htc) and Fig.[2tb)). 
Transmission rate in the new {y n } is given in Fig. (3Jb), 
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FIG. 3: The transmission rate T(e) of a 30-site chain coupled 
with metal electrodes with a voltage bias Vb = 0.48 (a), 0.56 
(b), and 0.72 (c). The coupling constant V = 0.32. The dashed 
lines indicate the Fermi levels of reservoirs. 



(a) 



(b) 



(c) 



(d) 



FIG. 4: Schematic energy-level diagram of a polymer chain 
for the formation of a soliton lattice, (a) A dimerized lattice 
with half-filled electronic band; (b) A dimerized lattice with 
one electronic excitation from a half-filled band; (c) One SS 
pair excited and (d) two SS pairs excited electronic band. 
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FIG. 5: Steady electronic current through a polymer chain 
of 30 sites (solid line). As a comparison, we also show the 
results for a uniform dimerized lattice (dashed line). 



creases to 0.72, the Fermi surface reaches the n = ±2 
levels of the polymer and another pair of solitons is ex- 
cited. When Vb increases again, soliton-latticc along with 
the corresponding soliton-lattice energy band is formed 
(Fig. E|d) and Fig.^c)). From the results we presented 
here, we conclude that the electronic transport through 
a polymer chain is accomplished completely through the 
soliton-lattice energy band. 



C. Nonequilibrium Electronic Current 

In Fig. [5j we show the nonequilibrium electronic cur- 
rent though a chain with 30 sites. It's clear that when 
bias voltage V> is below a critical value about 0.6, the 
current / is zero, which indicates a gap about 0.6 width 
existing in the energy band under these bias voltages. 
And the step-like feature in current-voltage curve is a re- 
sult of the finite size effect and can be explained by the 
gap shifted to match the Fermi levels of reservoirs, as 
shown in Fig. [3] As a comparison, we have also shown 
the current through a chain of a "uniform dimerized" , 
where the staggered lattice order parameter is taken to 
be independent of the site, i.e. y n = y for all sites n. 
It can be seen that the step-like feature disappears for 
a uniform dimerized lattice. The reason is clearly seen 
from the transmission rate T(e) of a uniform dimerized 
chain, shown in Fig. |6l where a band gap reduction by 
increasing the bias voltage is clear seen while no band 
shift be seen. 



D. Nonequilibrium Peierls Transition (NEPT) 



from which you can see one peak is moved down due 
to the SS pair formation for the energy levels of the soli- 
tons located inside the gap. It can be seen that the gap is 
moved up and matched the Fermi levels of reservoirs. Ex- 
cept the gap, other vales in the transmission rate shown 
in Fig. [3] are caused by the finite size of the chain and 
should disappear at thermodynamic limit. When Vb in- 



In a one-dimensional system, the lattice is not sta- 
ble against a 2kp lattice distortion, which is known as 
the Peierls transition.— For example, the ground state 
of a half-filled electronic system, such as polyacetylene, 
should be with a dimerized lattice at the temperature 
below a critical temperature. Ajisaka, et a/.,— consid- 
ered a sandwich structure that is an open chain coupled 
with two metal electrodes, in which the chain is treated 
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FIG. 6: The transmission rate T(e) through a uniform dimer- 
ized chain of 30 sites for the bias voltage Vb = 0.6 (a) and 0.7 
(b). 



by the Takayama-Lin-Liu-Maki (TLM) modclj±£ a con- 
tinuum version of the SSH model while the metal elec- 
trodes are described by reservoirs of free electrons. For 
a "uniform dimerized" lattice, they showed analytically 
that when the whole system is standing in noncquilib- 
rium steady states, a NEPT between ordered and nor- 
mal phases is found induced by the bias voltage Vb- The 
lattice is dimerized at the bias voltage below a critical 
value; otherwise the lattice will be uniform, i.e. y = 0. 

Our analytical derivation for a uniform dimerized SSH 
chain at thermodynamic limit is given in the appendix 
C. No different result is found between the continucm 
TLM model and the discrete SSH model for the NESS. 
And numerical results for a lattice with finite sites are 
shown in Fig. [2a) . A jump of the order parameter y 
ocuurs at bias voltage about 0.6, which is the critical 
value that electrons from the reservoirs can flow through 
the polymer chain. As a comparison, we also show the 
result at thermodynamic limit, which seems to us the 
jump as a sign of the NEPT. 

In Fig. [7Jb), we show results for a lattice allowed to 
relax. Here we introduce a parameter that is defined 
as 



f 



N ■ 



iV-l 



(24) 
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be independent of n. 




Energy Level Index n 

FIG. 7: Degree of lattice distortion y of a polymer chain with 
20, 30, 40 and 50 sites under different bias voltages Vb and 
(a) the lattice is restricted to be uniform and (b) is allowed to 
relax, (c) Degree of lattice distortion of an isolated polymer 
chain (30 sites) with different occupation n, which indicates 
the levels between —n and +n being singly occupied, as is 
shown in Fig. [4] In (a), y of a polymer chain at thermo- 
dynamic limit is shown. Stable solution is represented by 
magenta line while unstable by orange line. 



isolated polymer chain with different occupation n, which 
indicates the levels between — n and +n being singly oc- 
cupied, as is shown in Fig. @] The consistent with the 
results shown in Fig. [7] (a) and (b) told us that the lat- 
tice configure at noncquilibrium steady state with a bias 
voltage should be in correspondence with that an iso- 
lated chain with electronic excitations. With this, we 
can conclude that the NEPT does not survive the lattice 
relaxation and the soliton excitations, since the electronic 
excitations in an isolated chain only cause soliton excita- 
tions and the formation of a soliton latticed 



to represent the degree of lattice distortion. It is clear 
from the figure that y could reduce to 0.03 at least, not 
zero, which is different from that for a uniform dimerized 
chain. At the same time, we still see a jump at the bias 
voltage 0.6. These are contradict signs for the NEPT. 

To clarify the NEPT for a chain with soliton excita- 
tion, we give the order parameter y, in Fig. [TJc), for an 



IV. CONCLUSIONS 

We have studied the steady current and the lattice 
configuration in a mctal/polymer/metal structure with 
the metal electrodes treated as infinite reservoirs and the 
polymer chain described by the discrete SSH model of 
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electron-lattice interactions. Since the reservoirs are infi- 
nite, the whole system can reach a nonrqilibrium steady 
state (NESS) with persistent particle flow. In this situa- 
tion, it is the scattering states to be complete eigenstates 
of the whole system and the scattering state operators 
satisfy the Fermi distribution. Using the scattering state 
operator method, we show that in NESS, the solitons 
could be excited. Nonequilibrium Peierls transition does 
not survive the soliton excitations and the soliton-latticc 
is formed at thermodynamic limit. And it is the soliton- 
lattice energy band, not the conduction or the valence 
band, that provides the channels for electronic transport. 
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Now we obtain 



c { n a \e) 



where 



-1 / UR) 



for n > 2 and c{ (e) = 1, and 



for n < N - 1 and c^(e) = I. 



(A4) 



(A5) 



(A6) 



Appendix A: Calculation of the functions F^\e 
Based on Eq.([T0]). pi]), and fig]), we denote 

Fi%{e) _ F { n R) {e) E^fe) 



(Al) 



for n = 1 , 2, . . . , N — I, then the self-energies Y,^ (e) and 
E^ (e) satisfy following iteration relations 



4 L) (e) 



e-Si fi) (e) 



, (A2) 



with 



Appendix B: Calculation of the current flow 

The current flow from the left reservoir to the polymer 
chain is given as 



(Bl) 



fcf7 



and similarly, the current from the site n of the chain to 
the site n + 1 is, 



In = ^„^Im(4 i0 .c n+ i !(7 ). (B2) 



S^^e) = -iT R , Z[ n> (e) = -iT L . (A3) Based on the Eq.flTD), (O, and (JI5JI, we have 



afc' 



-/a(eafe') 



= . 9Lfc pf^fcJ/iCeLfc) + *r L |F 1 (L) (e ifc )| 2 /L(£L fe ) + ir fl |F 1 (ii) (e ifc )| 2 /i?(eLfe)J , (B3) 
where we have used the definition of DOS of the reservoirs (|I2[) in the large bandwidth limit. Then we have 

^ = f v/ *{-M^ (L) W}/L(e)-r i |F 1 (L) ( £ )| 2 / i ( e )-r J? | J F 1 (i?) ( £ )| 2 /«( e )}, (B4) 

where we have used £) fc 9lkf( € ak) = J def(e) J2 k dlk S ( e ~ e afc) = ( T a,/n) J def(e) for any function /(e). Furthermore, 
we have Im 



F^ L \e) = -T L \F[ L \e)\ 2 - tilm F 1 (i) *(e)F 2 (i) (e) from Eq.flS) and {KQ, so that 
fc = |^ / de{tam[F^>(e)F^(e)]f L (e)-T R \F^(e)\ 2 f R (e)}. 



(B5) 
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Now we see that the above expression goes to Eq. (f2"0l) if 
we can show 



film 



F^(e)F^ L \e)] = T R \F^ R) (e)\ 2 . 



(B6) 



For that, from the Eq. (|Al[) . we have 



f„Im 



F^>(e)F^ 1 {e)\=-JmE^(e)\Fi L \e)f, (B7) 

together with Eq. (|A2|) . we have 

ImSW(e)|Ff)(e)| 2 = ImS^e)^^ 2 , (B8) 
then we get 



LJm 



Fi L) *(e)F r {l\(e 



t„_ilm 



Fi L _\*(e)Fi L \e) 



which leads directly to 



(B9) 



Fi L) *(e)Fi L + \(e) 



t N _ilm 



4 L i*(e)4 L) ( £ ) 



The equality 



(BIO) 



(Bll) 



is easy to be verified directly from the expressions of 
functions Fn a \e) obtained from Eq. (|A4[) . which does not 
require the central inversion symmetry of the coupling 
with the reservoirs and the lattice configuration. Then 
we proved Eq. (lB6p . 
Similarly, we have 

(ct i(T c„ +1>CT > = Y,\9ak\ 2 F 7 [ a >(e ak )F^ 1 (e ak )f a (e ak ) 

ak 

a J 

(B12) 



so that the current of the Eq. (|B2[) becomes 
2e 2 



^r a y de t„Im [F^*(e)F^ 1 (e)\ f a (e). 

(B13) 



Doing the same as Eq. (|B9[) . we have 



f„,Im 



film 



Ff)*(e)F 2 ^)(e) 



-rx|Fr ; (e 



(B14) 



together with the Eq. (|B10[) and (|B11[) , we show that /„ 
is independent of site index n, and finally I n = II = I, 
which meets the requirement of a NESS. 



Appendix C: The case of a uniform dimerized lattice 
at thermodynamic limit 

From Appendix B, we rewrite the transmission rate in 
Eq. tpTj) as follows 



T( e ) = -4r i |FW( £ )| 2 ImEW( e ). (CI) 
From Ea.(|A"4)l together with Ea.(|AT]l and (|A2|) , we have 
1 ImS[ fl) (e) 



\Fk L \e)[< 



rL| £ -sW( e )-S^( £ )l2 : 



then we have 



T(e) 



4ImS^ ) (e)ImSr j (e) 



(R)t 



(R), 



(C2) 



(C3) 



which has been shown to be a quantity independent of 
index n in the NESS. 

The bond-charge density p n , n +i given in Eq.(f2U]) can 
be written as 

P^n+^E/^W/^ 6 )' ( C4 ) 



where 



{L) 4 Regg(e)IrnEg)(e) 



£ -C'(£)-£[ fl, (e)P 



and 



^ } (e) 



4 Rgggi 



(C6) 



For a polymer chain of a uniform dimerized lattice at 
the thermodynamic limit, y n = y, the hopping constant 
t n = 1 + (— l) n y (hereafter, to(= 1) as the unit of en- 
ergies), and the self-energies S„(e) for site n being far 
away from both of the ends should satisfy the following 
self-consistant equations 



4 L) (e) 



t: 



, Sn+i(e) — 



rt+l 



e - S„(e) 



(C7) 



where we have used £„ 



(i) 



J n+2' 



Then, we have 



ImS„(e) = for |e| > 2 or |e| < 2\y\, while for 2\y\ < 
lei <2, 



^) (£) = £ a + (-l)"4 g _.V(4-y-V) (C8) 



2e 



2d 



Similarly, sl fl '(e) is obtained to be that of (e) re- 
placed y by —y. From the self-energies we obtain the 
transmission rate T(e) = for |e| > 2 or |e| < 2\y\ and 



9 



T(e) - 1 for 2\y\ < \e\ < 2. And p^\e) = p ( n R) {e) 
p„(e) as given below 



e 2 + (-l) n 4j/ 



|e|7rt„ ^(4- £ 2)( e 2_ 4y 2 ) 



(C9) 



for 2 1 y | < |e| < 2 and vanishes otherwise. 

The lattice dimerization y is determined by f(y) = 
with the function f(y) is given as 



/(») = ^-^E(- 1 ) n ^.M-i 



(CIO) 



where we have used the Hcllmann-Fcynman theorem. In 
the thermodynamic limit, the contribution of the sites 
close to the two ends, where the bond-charge density 
might be affected by the reservoirs, is negligible. Then 
we have 

f(y) = -T - \{P2n.2n+l ~ p2n+l,2n+2) (Cll) 

with 2n being the site far away from the two ends so that 
the right hand of the above equation will be independent 
of index n. It's clear that the lattice dimerization y will 



be independent of the couplings with reservoirs since the 
dominant bonds in thermodynamic limit are those far 
away from the two ends. 

For a symmetric applied bias voltage fir, = — A*r = 
Vb/2 with 2 > Vb/2 > 2\y\, the self-consistant equation 
for y becomes 



Ay 



i - v 2 J- 



-V b /2 



de 



4y 2 



(C12) 



The integration from —Vb/2 to Vb/2 in the above equa- 
tion for /i(e) has been cancelled out due to the anti- 
symmetry of the integrand. For the bias voltage Vb < 
4\y\, the integration in Eq.(C12) will be up to — 2\y\ and 
the self-consistant equation for y will return back to the 
one for an isolated chain, independent of applied bias 
voltage Vb- 

The lattice order parameter y under different bias volt- 
age Vb can be obtained from Eq. (|C12[) . shown in Fig.lTJa), 
from which a NEPT can be seen as that of the TLM 
model obtained by Ajisaka, et al.^ No different is found 
between the continuem TLM model and the discrete SSH 
model for the NESS. 



* Email: cqw@fudan.edu.cn 
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